Testing Different Formulations of MOND Using LISA Pathfinder 
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Previously it has been shown that Modified Newtonian Dynamics (MOND) can be tested near 
the saddle points of the Newtonian gravitational potential using the forthcoming LISA Pathfinder 
(LPF) mission. These analysis focused on one particular formulation of the MONDian theory, 
which centered around a particular modified Poisson equation. We show that, in addition to the 
well known AQUAL formulation, another possibility exists in the form of a driven Poisson equation 
for the MONDian field. We consider similar quantitative and qualitative analysis for this formulation 
and also investigate typical Signal to Noise Ratios (SNR) resulting from these theories for a LPF 
test. We demonstrate that for a typical 50 km saddle flyby, the SNR is amplified by 25% between 
the two formulations. We also suggest that the SNR remains enhanced for impact parameters of 
1000km or larger. Finally we show that constraints from a negative result remain as good, but no 
better between different formulations. 

PACS numbers: 04.50.Kd, 04.80.Cc 
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I. INTRODUCTION 

In the past few decades, attempts to reconcile anoma- 
lous observations of galactic dynamics without resorting 
to dark matter have been one of the driving forces behind 
modified gravity theories. These theories have centred 
around using a prescription for a modified force law 
which suggests the usual Newtonian dynamics at 'large' 
accelerations and modified dynamics at 'small' accelera- 
tions (large and small here are relative to the Milgrom 
characteristic acceleration ao ~ 10~ 10 ms~ 2 ). In doing 
so, the observed flattening of galaxy rotation curves and 
the empirical Tully-Fisher relation can both be satisfied. 

Recent work [2 4] on an experimental observation of 
MONDian like theories suggests that anomalously large 
tidal stresses should be present around gravitational sad- 
dle points. The proposed LISA Pathfinder (LPF) mis- 
sion Q should be capable of probing acceleration scales 
upto ao, this is of special importance because here we 
would be probing the onset of MONDian behaviour, 
rather than the acceleration scales dominated by MOND. 
The purpose of this short paper is extend previous anal- 
ysis done of the MOND saddle point science case, based 
upon a scenario where a mission extension is granted. 
The extension would involve redirecting the spacecraft 
from Lagrange point LI to a saddle of the Earth-Moon- 
Sun system ,(| once its nominal mission at LI is com- 
pleted. Previous analysis have only considered only non 
rclativistic limits resulting from theories such as TeVeS 
or ones which are AQUAL- like, here we will examine how 
these predictions vary using instead a different formula- 
tion of MOND. 

The structure of this paper is as follows, first we re- 
view the non-relativistic limits of the various modified 
gravity theories in the literature. In section III we solve 



the Poisson equation for this theory numerically around 
the Earth-Sun saddle point. We consider typical signal 
to noise ratios for detecting MONDian effects with LISA 
Pathfinder in section V, our main results being sum- 
marised in figures |3] and |4] Some more general concerns 
regarding these type II theories are considered in the last 
section and we will conclude with some thoughts on using 
all of these results constraining modified gravity theories. 
We also present some analytical solutions and detail our 
numerical methods in appendices. 



II. A REVIEW OF MONDIAN THEORIES 

The literature contains a whole plethora of relativistic 
MONDian theories, but we find their relative differences 
(and somewhat complexity) arise from the requirement 
that they explain relativistic phenomena (such as lensing 
and structure formation) without resorting to dark mat- 
ter. In the non-relativistic regime however, all of them 
reduce to just to only 3 types of non-relativistic limits, 
which we will label type I, II and III. 

• Type I In these theories, the non-relativistic dy- 
namics results from the joint action of the usual 
Newtonian potential V$at (derived from the met- 
ric via goo « — (1 + 2V$at)) and a fifth force field, 
4>, responsible for MONDian effects. The total po- 
tential acting on non-relativistic particles is their 
sum: 



$ = $at + (f) 



(1) 



Whilst the Newtonian potential satisfies the usual 
Poisson equation: 

V 2 $ w = AnGp (2) 

the field <j) is ruled by a non-linear Poisson equation: 
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V • (mO)V<£) = kGp 



(3) 



2 



where, for convenience, we pick the argument of the 
free function function \i as: 



k |V0| 
Air clq 



(4) 



where k is a dimensionless constant and ao is the 
usual MOND acceleration. In general, we require 
that [i — > 1 for z ^> 1 and ^ ~ z for z < 1. 

• Type II In these theories we also have $ = &N+<f), 
but now the field is ruled by a driven linear Pois- 
son equation, whose source depends on the New- 
tonian potential. In order to facilitate comparison 
with Type I theories, we write the equation for cj) 
in these theories as: 
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V ■ (v(v)V<Z> N ) 



(5) 



where the argument of free function v is given by 



|V$ 



N 



c 



(6) 



and we require that v — > 1 when v 3> 1 and 
2/ ~ l/v 7 ^ for w < 1. It is the analysis of the 
anomalous tidal stresses around the Earth-Sun sad- 
dle in this non-relativistic limit that is the subject 
of this paper. 

Type III This was the original non-relativistic 
MONDian proposal, derived from a non-relativistic 
action principle (the so-called AQUAL 7]). Cru- 
cially, here non-relativistic particles are sensitive to 
a single field $ which satisfies a non-linear Poisson 
equation: 



V • (AO)V$) = AirGp 



(7) 



Again, \i is a free function with a suitably chosen 
argument: 



|V$| 

x = 

a 

so that fi — > 1 for i» 1 and jl ~ x for x -C 1. 



(8) 



Virtually all relativistic MONDian theories proposed 
in the literature fall into these categories. Both 
TeVeS, Bekenstein's relativistic MONDian theory [8j and 
Sanders stratified theory 19 1 have a type I limit. Mil- 
grom's bimetric theory [JiJ, LUj can have either a type I 
or type II limit, whilst Einstein- Aether theories [Jja, [U| 
appear to be to the only theories with a non-relativistic 
limit of type III. 

We stress that it is important that we consider for- 
mulations wedded to fully relativistic theories (rather 
than just a scalar Lagrangian theory like AQUAL), even 
though we are not directly testing these full theories 
themselves. Their cosmologies however, do provide the 
necessary constraints on the gravitational constant G, 
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FIG. 1: A comparison between the numerical and 
analytical results for components of g = — V0 in the 
two-body Earth Sun case. 



the bare value of G appears in the friedman equations 
and so from BBN constraints can be fixed to corrections 
of order ~ 0.01 or smaller. As we will see, this issue is 
not a trivial one in modified gravity theories. Deriving 
the full field equations and taking the weak field limit as 
necessary allows us to find the modified Poisson equa- 
tions for these theories consistently. A more thorough 
investigation of the weak field limit in these theoreis is 
left for future work [3] 

It has been argued that in some relativistic formula- 
tions of type II theories the bare G, that appearing in cos- 
mology and the total G ruling the non-relativistic equa- 
tion are the same. We call such theories type IIA, and 
for them v — > in the Newtonian limit. Otherwise we 
call then type IIB, with a G renormalization, and with 
v^t 1. 



III. NUMERICAL SOLUTIONS 

To make further progress with the saddle point prob- 
lem, we must move to using numerical solutions. We use 
a modification of the code used by Bevis et al (see Ap- 
pendix [B] for details), which solved the modified Poisson 
equation over a lattice around the saddle. Non-uniform 
coordinates are employed to aid resolution in the deep 
MONDian regime. Due to the linear nature of equation 
([5]) , a variety of computation techniques can be used (no- 
tably FFT). Here however we have the advantage of using 
our code, whos solar system credentials have been exam- 
ined previously Q . Solving the sourced Poisson equation 
on the lattice to find SF — — V0, we find similar results 
to the Type I theories. We picked a v function of the 
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form 



1 + - 

v 



1/4 



(9) 



since it has similar properties to the /i function picked for 
type I theories, however it is simple to incorporate other 
functions in our code. 

As before, we considered the Earth-Sun saddle (al- 
though three body extensions can be easily considered 
too), using a 257 3 lattice of physical size 10000 km with 
central resolution ps 2.6 km. As figure [T] shows, these pro- 
duce a good match to analytical solutions (within their 
respective domains) and provide the appropriate inter- 
polation in the intermediate region around the bubble 
boundary. 



IV. TIDAL STRESSES AND SIGNAL TO NOISE 
RATIOS 

Whilst predictions of MONDian forces are useful, LPF 
is sensitive to the tidal stresses. In addition to the is- 
sues regarding close flybys for the spacecraft (see [(| for 
a further account), two additional ones need to be taken 
into consideration here. One is that only transverse tidal 
stresses can be measured, due to practical issues regard- 
ing the spacecraft. The other is that (f> field produces 
both a MONDian effective field and a rescaled Newto- 
nian component, which both must be taken into account 
when calculating the tidal stresses, so we use the defini- 
tion: 
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dxidxj Att dxidxj 



(10) 



As stressed before, it is important that we know both field 
components to the same degree of accuracy or systematic 
errors can crop up due to the imperfect subtraction of the 
Newtonian field (see [H, section IV C). 

We then move on to consider typical signal to noise 
ratios (SNR) produced by LPF in this formulation of 
MOND, using the techniques of noise matched filtering 
from gravitational wave searches fl5| . The basic idea 
laid (as laid out in 4]) is to correlate a time series x(t) 
with an optimized template designed to provide maximal 
SNR, given the signal shape h[t) and the noise properties 
of the instrument. Whilst the specifics of the instrument 
noise won't be known properly until the satellite is in 
situ, there exist nominal requirements that it must meet, 
as well as best estimates for the noise signal waveform. 

For the setup we described here, we have h(t) = 
S vy (vt,b,0), where v is the velocity of the spacecraft 
(taken to be 1.5 km" 1 throughout this work) and t = 
corresponds to the point of closest saddle approach. In a 
more general setup, for an approximately constant veloc- 
ity v, a closest approach vector b, and with the masses 
aligned along unit vector n, we have 
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FIG. 2: ASD plot of the MONDian and Newtonian 
signals (multiplied by k/4tt), as compared to the basic 
noise ASD, for both the Type I and Type II theories, 
for a y = 50km fly-by, with v = 1.5 ms" 1 . Using this 
data, we find the SNR is amplified from 28 to 35. 



The maximal SNR, realized by the optimal template, is 
found to be: 



Popt 



4f 



Kf) 



Shif) 



1/2 



(12) 



h{i) = n i n j S ij {h + vt) 



(11) 



where h(f) is the fourier transformed signal waveform, 
computed from the tidal stresses and Sh{f) is the noise 
waveform. We plot these in figure [2] (in terms of am- 
plitude spectral density (ASD), ie square root of power 
spectrum) to demonstrate that there is ample signal com- 
pared to the noise for LPF to detect MONDian behav- 
ior for a reasonable fly by. The headline figure is for a 
b = 50km flyby, with baseline noise 1.5xl0~ 13 s~ 2 /\/Hz, 
the SNR observed should be « 35, a distrinct improve- 
ment on the figure of f=s 28 using the Type I theories. 

As before, we can produce contours for typical SNRs 
produced at various impact parameters and baseline 
noises, see figure [3] In doing so, we suggest that larger 
SNRs (in comparison to type I theories) will in general 
result, for b's in excess of 1000km, as illustrated in figure 
21 The most striking thing to take away from this plot 
is that even as far out as 1000km, there is still a marked 
difference between SNRs of the two theories. 



V. FEATURES OF TYPE II THEORIES 

Given the amplified MONDian signal that we showed 
was present in these theories, it is tempting to suggest 
that a null result would represent an even tighter con- 
straint on these theories than compared to type I. We 
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FIG. 3: Signal to noise ratio (SNR) contours for various 
impact parameters and baseline ASD noise. We set the 
v = 1.5 kms" 1 . Calamitous assumptions would still lead 
to SNR in excess of 5. More optimistic ones (b around 
50km or less, noise half way up the scale) would lead to 
SNRs easily around 55. 
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FIG. 4: Signal to noise ratio (SNR) contour lines for 
various impact parameters and baseline ASD noise 
waveforms comparing between type I and II theories. 
The solid lines are the typical SNR to be obtained in the 
type II theories and the dashed lines to their immediate 
left the corresponding type I line. We show for impact 
parameters up to b = 1000 km to demonstrate type II 
always substantially beats type I on SNR. 



should however be clear to distinguish between a stronger 
MONDian signal and the relative size of the Deep MON- 
Dian region tq. In these theories, with 'simple' single 
power law v functions, the bubble size ro remains un- 
changed, this is obvious (eg from analytical solutions, see 
appendix [A]) . The tidal stresses however are amplified at 
larger frequencies (see figure [2]) , the integrated effect of 
which is a larger SNR. This however is not the same thing 
as a theory with a larger MOND bubble, which would 
boost SNRs at all frequencies (equivalent to making a 
closer fiyby). 

A. One or Two Field Models 

We observe that §5§ can be rewritten as a manifestly 
one field model in terms of the gravitational potential $ 
and an (auxilliary) Newtonian Potential i'jv as 

V 2 $ = V ■ (f>(w)V$iv) (13) 

where 

Kw) = 1 + |^(«) ( 14 ) 

and w — | V$ at | /do- An observation we can make in this 
case is, that for large Fn, u — > 1, making 

G Ren = (l + A) G (15) 

whereas if we picked a manifestly one field model (HA), 
there would be no G renormalisation and we would have 
v — >■ in the large Fjq limit. This has the entertaining 
consequence of being equivalent to a model with /x — > oo, 
a result we discuss in the next section. 

B. Divergent functions - a discussion 

In Type I theories, we have the case where 

V-(mF ) = Ay.Fjv (16) 

which under the assumption of spherical symmetry, we 
can ignore curl terms and reduce to 

/iF^AFjv (17) 

We find there are broadly two classes of [i functions, those 
which converge to some constant and those which diverge 
in the large acceleration limit. 

Lets consider the sub class ofu's parameterised by a, 
with a = 1 first suggested in [l6| as being a good fit to 
galactic data: 



(18) 
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where for brevity k' = k/A-K. Here it appears that for 
z — > oo, /i — > — ~, however it turns out if we investigate 
the behaviour of z properly, by substituting into (|17l) 

(19) 



1 fc' 



and then solving for z (picking solutions for z > 0) 
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4(fe 



A2 



and so as F, 



N 



T7 a F W 



(20) 



(21) 



i.e. z goes to a constant as /x diverges and we never reach 
the other branch of solutions for /i(z > ~). This also 
means that since the MONDian field approaches ~ a 
asymptotically, we are always in the region where 'deep' 
MOND effects are present (see section VIA of [4[). We 
can however use (12011 to write this asav function 



a 

2k' 



4(fc') 



(22) 



which evidently has v — > + in the large Fjq limit. This 
has the predictable conclusion that with these functions, 
Grsh = G. It is however exactly these functions which 
are required in a one field model. 

We can schematically write a function with the same 
behaviour as 



(23) 



1 



ie the usual v — > m the MONDian regime but 

moving to a different power law for larger accelerations 
v — > I /v. This notion of switching power laws at different 
acceleration scales is the subject of the next section. 



C. Designer v functions 

One requirement of an experimental test of MONDian 
theories is to either conclusively show there is an extra 
gravitational force present in the low acceleration regime 
or in the case of a null result, provide clear constraints 
for theories. It is reasonable now to suggest that fly by 
impact parameters that can probe the MOND bubble (ie 
r < 380km) will be possible. With this in mind, lets 
consider probing that regime and the effect of seeing no 
signal above the Newtonian background and the noise. 
At best one, can then posit SNR of order 1 and the size of 
the bubble is smaller than originally thought. We would 
then have to design a transition function as such to escape 
the net of LPFs capabilities in this way. We considered 



this approach with the /z functions of type I (see section 
VI of [J] ) and here we apply it similarly to the v functions 
of type II. 

We start by fixing the asymptotica of our problem, 
which are v « 1/ \fv in the astrophysical regime (ie 
Fn < o-o) an d v ~ 1 in the quasi- Newtonian regime 
as such to be washed out far away from the saddle by 
superior size of the Newtonian potential. In the inter- 
mediate regime (which we are probing with LPF), we 
can suggest a different power law, as such making this a 
two scale model, with the power n and ao now the two 
independent features of these theories to be constrained. 
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(24) 

trig ( 25 ) 



trig 



(26) 



where the point when non-Newtonian behaviour in 
triggered can be interchangeably pinpointed by: 



..trig _ 



trig 



trig 

1 N ~ 



f-V 



0() 



(— 



(47J 



(27) 
(28) 
(29) 



We still have that when < 0,0 the field <fi dominates 
$jv as per our requirements, but now the intermediate 
region, where <f> hasn't yet dominated but is already non- 
Newtonian, is in a narrower band of accelerations ao < 



trig 
N 



As a result, the MOND bubble shrinks in 



this model according to 



ro 



380 

V47T 



km 



(30) 



Making a model independent statement on which power 
would be required for a SNR of order 1 remains beyond 
our reach, simply because a null result can only lead us 
to conclude we are probing the region b 3> ro(n). Equiva- 
lently this statement means we are probing the transient 
from truely (albeit rescaled) Newtonian behaviour v ~ 1 
and power law behaviour v as v n . This however can vary 
wildly even between 'simple' v functions of one power. 

We can however perform an order of magnitude argu- 
ment, assuming a designer function, say: 



v(v) = 1 + 



trig 



(31) 



1 Where again we notice that the smallest the bubble can shrink 
to in this model is ro ~ 1km. 
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FIG. 5: Contours of the power n needed to obtain SNR 
= 1 for different impact parameters and noise levels, 
upto b — 600 km For n > 2, we find our functions 
become increasingly contrived in order to explain a null 
result. This suggests that a negative result at ro/2 or 
closer would need to be explained away differently. 



We can write 

F <t> = ^uF N (32) 

47T 

since it is reasonable to ignore divergence free sources 
here (y « 1 at these scales). Breaking this up into 

F = °F + SF^ (33) 

with rescaled Newtonian component "F^ = j^F^r- Sub- 
stituting in and solving gives: 

„,_*<,,_,)„„_ (34) 

from which the tidal stresses can be inferred. We notice 
that in this instance, this equation is exact since v = 
v(Fn), whereas for type I theories, this process was done 
perturbatively. 

The results of using this model are presented in figure 
[5] (kindly borrowed from [J]), indicating the value of n 
needed for a given b and noise level in order for a SNR 
of one to be obtained (and so a negative result be ac- 
ceptable). Once we start probing deep into the MOND 
bubble, a rather 'fine tuned' designer v becomes neces- 
sary to accommodate a negative result. The naive ex- 
pectation that these null results are more constraining 
on type II theories is this a misnomer, because we are 
not probing the deep MONDian regime here, merely the 
transient from quasi-Newtonian behaviour. In this way, 
a null result would make differentiating between theories 
inconclusive. 



VI. CONCLUSIONS 

In this paper we solved the non-relativistic limit of 
a different MONDian theory formulation (here classed 
Type II) for the saddle points in the two-body Earth- 
Sun system. We probed the additional MONDian effects 
and forecasted for LPF flyby probable detection possi- 
bilities. We suggest that the chances of detection of a 
MONDian signal are enhanced compared to Type I the- 
ories (eg those resulting from TeVeS), with typical SNRs 
for a 50km flyby rising from 28 to 35. Another key result 
is that the boosted SNRs remain higher even as far out as 
b « 1000km or more, enhancing observational prospects 
even for catastrophic impact parameters. We also sug- 
gest that in the event no anomalous tidal stresses are 
observed, the experiment can put constraints on these 
theories, however they remain no better than in Type 
I, due to the nature of the regime we would be prob- 
ing. Further work however remains to be done into ways 
to differentiate between MONDian formulations. One 
significant difference here is that there are no curl type 
forces present, due to the linear nature of the governing 
equation. Whilst a null result would place constraints 
on these theories, a positive one would do little to guide 
us towards a particular formulation or even form fi or v 
function. We are currently working on techniques to con- 
strain or perhaps extract some features of the function 
from potential data obtained. 
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Appendix A: Analytical Solutions 

In this work we will follow the conventions and no- 
tations used previously. Our experimental observations 
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to ours, ours is a two field model. They on the other hand use 
a one field model, where the MONDian field is an auxiliary one 
and the observable is the gravitational potential <3>. Our tech- 
niques are kept similar to previous work to allow more direct 
comparison between the two theories. The reason that the two 
approaches are so complementary for this test, is simply that we 
are free of sources at the saddle, so the resulting equations are 
the same, this result would not be in general true. 
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will be centred around the saddle point so it makes sense 
that our coordinate system is too. For solving the Pois- 
son equation about the saddle, we pick a radial system 
of coordinates, however when we consider spacecraft tra- 
jectories we revert to Cartesian coordinates. 

Although we have freedom in picking the form of the 
v function, the one used throughout this work was 



" = 1 1 + 4 

V 



1/4 



(Al) 



This was picked as so to easily compare with previous 
results from type I theories. 

About the saddle point, we have the Newtonian field 



where 



V$at = g N = ArN 



N(ip) = N r e r + N^ 
N r = i[l + 3cos(2V0] 



■sin(2^). 



and hence 



v = —N 



(A2) 

(A3) 
(A4) 

(A5) 
(A6) 



where N = |N| ~ 0(1). 

Replacing the Newtonian force in ([5]) with (|A2|) , using 
the fact that at the saddle V 2 $at =0 and expanding out 
gives our expression in this linear regime: 



<f) = Y (t 



1/2 



1/4 



r o r ) V r o + (rN) 2 

N r . Nj, dN 
N 1 / 2 + N 3 / 2 ~£hp 



(A7) 



The problem here is akin to electrostatics, where we must 
solve the sourced Poisson equation subject to the ap- 
propriate boundary conditions, which here are that SF^p 
vanishes (and 5F r equate) at ip = and 7T, such that we 
avoid a jump in the field at -0 = ir/2 Q. 



1. 'Deep MONDian' 

For r -C r , it's clear (|A7I) reduces to: 

72 _ a /4tt\ / 1 \ 1/2 / N r N^, dN 



2 \ k J \rr Q J V^ 1/2 N3/2 d i> 



(A8) 



where the angular parts of the leading order term neatly 
reduce to 



7 + 9 cos 2tfj 
(2(5 + 3 cos 2^)f/ 4 



S«0 



(A9) 



and the higher powers are supressed since »<1. 
We use an Ansatz for the MONDian field 

4> = Ci r a F(il>) 



(A10) 



where C\ is some constant to be fixed from the source 
term above. This gives rise to a sourced 2nd order ODE: 

r a - 2 (a(a + l)F + cot(^) F' + F") = r~ 1/2 g(tp) (All) 

where ' = 8/ dtp and 

47T an 



(A12) 



We find that the solutions of the homogenous (un- 
sourced) equation are Legendre Polynomials of order a, 
with the form of (|A11|) suggesting a — 3/2. The full in- 
homogeous (sourced) solution is then given by the vari- 
ation of parameters method and the solutions given as a 
Fourier series expansion in even components of tp: 

F{i))K - 0.0236032- 0.188606 cos(2^) (A13) 
+ 0.010834267 cos(4^) + ... 

In addition, at the saddle we have region where |gjy| = 
and so we need to consider solutions to the Laplace 
equation 







(A14) 



which subject to smoothness and continuity conditions 
being satisfied and regularity at the origin, can be written 
in general by 



4>L = a-o 



in 



J2 A e r t P 2 i{cos%p) 



where P2i(cosip) are Legendre polynomials, 



(if 



(A15) 



(A16) 



and ai are dimensionless constants to be found by match- 
ing solutions at the intermediate MONDian regime. Our 
normalisation is picked to be of the same form as the 
sourced solutions, so that V(j> has units of acceleration. 
It is important to consider that we need only expand a 
few terms out form the Laplace contribution when com- 
paring against data, since the region of validity of these 
solutions is small. 



2. 'Quasi-Newtonian' 



For r 3> ro, it's clear that (IA7|) reduces to: 



ao -r 



Att\ ro 1 



k J r 2 2N 2 



aiVAv, 
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(A17) 
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where the leading order angular term simplifies to be- 
come 



2(7 + 9 cos 2ip) 
(5 + 3cos2^) 2 



h(iP) 



(A18) 



and again higher powers are suppressed since v» 1. 
Our Ansatz for the MONDian field here is 



h = C 1 H 2 (t/>)+C 2 \n — 



(A19) 



since if used a similar ansataz as before, we wouldn't 
be able to satisfy smoothness and continuity on [0, 7r] 
(we can however use them for higher order terms in the 
expansion of v). Computing the Laplacian gives 

r,2 > Ci 1 d , ■ C*2 /4tt\ r 0l , , N 

V 2 = sinff^ + -f r = ao — V 

r z sm ip oip r z \ k J r 



allowing us to set 



Ci = C 2 r = 

Integrating out once then gives 
dH 2 



4tt 



sin ip 



dtp 



(h- 1) sini/i dip + A 



(A20) 
(A21) 

(A22) 



Again requiring that we satisfy our boundary conditions 
means we find 



.4 



3 7T \ 

2 + 37!j 



(A23) 



We can then solve (|A22[) for H 2 and expand as a fourier 
series in ip 

H(ip)& - 0.22921 + 0.287573 cos 2ip (A24) 
- 0.116339 cos 4V> + ... 



Further terms in the expansion of (|AT7|) result m a 
series expansion for the MONDian field <p 



-in 



ao J]] Cn 

n=2 



2-2n 



H n {^) (A25) 



where H n (ip) satisfies the sourced ODE 

n(n + l)H n + cot{ip)H' n + H„ = K 
and h n is given by 

2 3 "/ 2 - 2 (7 + 9 cos 2^) 



h n {lp) 



(5 + 3cos2V>)™ /2+1 



(A26) 



(A27) 



In addition to these solutions, we also always have the 
background rescaled Newtonian contribution 

$Ar = ^f^-(l + 3cos2^) (A28) 
8 \ k J r 



In the limit r oo, we recover this rescaled Newtonian 
limit. 3 



Appendix B: Details of Numerical Method 

Here we describe our numerical algorithm. We con- 
sider an adaptation of the the numerical set up used 
previously Q, representing g = — Vcp on the sites of a 
non-uniform lattice around the saddle point. We frame 
the system as a pair of vector equations 



v .g = Jiv-tKffV] 

V x g = 



Our initial conditions are at each site set as 



N 



g0= 4^ g 



we then compute the discrete divergence 



x A3 



N 



where we will use the compact notation 

^\ ^ <Y* 3 IY*3 

- ~~ X x-j 
M - J 



and have discrete source term 



D 



N 



E 

3 



MA 



(Bl) 
(B2) 



(B3) 



(B4) 



(B5) 
(B6) 



(B7) 



such that at each site we locally solve (|B1[) whilst ensur- 
ing (|B2|) is satisfied globally. As before, we can consider 
changes in the field g J — > g J + Sg^ that keep the discrete 



(B8) 
(B9) 



curl condition (|B2[) satisfied, where 

»* A 3 + 



These changes are then reflected when we compute the 
change in the discrete divergence 5D X 



3 The choice of background terms is not unique, but rather depen- 
dent on the form of v. If we had picked a v with both odd and 
even powers in its expansion (i.e. v ss 1 + cei/v + ... ), there 
would also be a background 4>i contribution since V0i would be 
an angular function only. However this terms contribution in the 
r — > oo limit is negligible. 
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Ensuring that the change in the divergence at each step 
now gives 

Ac + 5D X = (BIO) 

meaning that the change in the MONDian field at each 
site is now given by 

« = -^t (B11) 

such that as we cycle through the lattice and the g 



field converges, the additional changes to g 3 x lessen. We 
achieve faster convergence using a successive over relax- 
ation method (SOR) as before, by scaling the field as 

SgL -+ (B12) 

where A is the over-relaxation parameter and is larger 
than unity. We begin with A = 1 and increase it once the 
field is settling down, since high values of A can initially 
result in the RMS value of \SD X \ increasing, whilst we 
arc looking for \5D X \ — > 0. 
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